% function constantsReal 

%% FOR FOCM side 

%stateNames=stanames; 
%parnames=coefnames; 
%parMode = param;
%Ynames = Znames;

[GG, RR, CONS, eu, SDX, ZZ, initss, ssvec,~, ssNames,stateNames,shockNamesStru]...
    =feval(funcmod,parMode,solveopt,addsol);
[KFStru,logL]=feval(addsol.funcKfilter,parMode,funcmod,Y,trainvec,solveopt,addsol);

%% FFR
posFFR=cellposition('FFR',stateNames);
paramFFR=CONS(posFFR);

%% Core PCE
paramInflation=parMode( cellposition('pi^{ss}',parnames) );
paramPCECore=CONS( cellposition('PCE Core',stateNames) );

%% Hours
paramHours=CONS( cellposition('Hours',stateNames) );

%% Real Average
paramGAM100=parMode(cellposition('gamma_{\\star_{100}}',parnames) );

%% GDP
paramGDPDef=CONS( cellposition('GDP deflator',stateNames) );
paramGDP=paramGAM100+paramInflation-paramGDPDef;

%% Consumption
paramCONS=paramGAM100+( paramInflation-CONS(cellposition('C deflator',stateNames)));

%% Investment
paramINV=paramGAM100+( paramInflation-CONS(cellposition('Inv deflator',stateNames)));



% findData=cell(3,2);
% findData(1,:)={'Nominal Per Capita GDP Growth','Inflation (GDP Deflator)'};
% findData(2,:)={'Nominal Per Capita Consumption Growth','Inflation (Consumption of ND & S)'};
% findData(3,:)={'Nominal Per Capita Investment Growth','Inflation (Durable Consumption & Private Investment)'};
% findStates={'GDP','Consumption','Investment'};
% findCONS=[paramGDP;paramCONS;paramINV];
% for ii=1:3;
%     Yreal=Y(:,cellposition(findData{ii,1},Ynames) )-Y(:,cellposition(findData{ii,2},Ynames) );
%     Ystate=KFStru.smoothSt(:,cellposition(findStates{ii},stateNames) )+findCONS(ii);
%     comparemat(Yreal,Ystate);
% end
    